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A Robust Distance Measurement and Dark Energy Constraints 
from the Spherically-Averaged Correlation Function of 
Sloan Digital Sky Survey Luminous Red Galaxies 



(N 

o 

(N 
X> 

IX, 
(N 



O 

u 

6 



> 

(N 
(N 
00 

00 

o 
o 



Chia-Hsun Chuang*, Yun Wang, and Maddumage Don P. Hemantha 

Homer L. Dodge Department of Physics & Astronomy, Univ. of Oklahoma, 440 W Brooks St., Norman, OK 73019, U.S.A. 



29 February 2012 



ABSTRACT 

We measure the effective distance to z — 0.35, IV (0.35) from the overall shape of the 
spherically-averaged two-point correlation function of the Sloan Digital Sky Survey (SDSS) 
Data Release 7 (DR7) luminous red galaxy (LRG) sample. We find IV (0.35) = 1428^3 
without assuming a dark energy model or a flat Universe. We find that the derived measure- 
ment of r s (zd) / Dy(0.3b) — 0.1143 ± 0.0030 (the ratio of the sound horizon at the drag 
epoch to the effective distance to z = 0.35) is more tightly constrained and more robust with 
respect to possible systematic effects. It is also nearly uncorrelated with fl m h 2 . 

Combining our results with the cosmic microwave background and supernova data, 



we obtain fit. = -0.0032 



+0.0074 
-0.0072 



and 



= -1.010 



-0.046 



(assuming a constant dark en- 



ergy equation of state). By scaling the spherically-averaged correlation function, we find 



the Hubble parameter _ff(0.35) 



83 



+ 13 
-15 



km s J Mpc 1 and the angular diameter distance 



Da (0.35) = 1089I|? Mpc. 

We use LasDamas SDSS mock catalogs to compute the covariance matrix of the corre- 
lation function, and investigate the use of lognormal catalogs as an alternative. We find that 
the input correlation function can be accurately recovered from lognormal catalogs, although 
they give larger errors on all scales (especially on small scales) compared to the mock catalogs 
derived from cosmological N-body simulations. 

Key words: cosmology: observations, distance scale, large-scale structure of Universe 



1 INTRODUCTION 

Galaxy redshift surveys provide a cosmological probe highly 
complementary to the cosmic microwave ba ckground ( CMB) 
(Penzias & Wils on! 19651) and supernovae (SNe) iRies s et al J 19981 ; 
IPerlmutter et all 1 999). Large-scale structure data from galaxy sur- 
veys can be analyzed using either the power spectrum analysis 
or the correlation function analysis. Although these two meth- 
ods are simple Fourier transforms of one another, the analysis 
processes are quite different and the results cannot be converted 
with Fourier transform directly because of the finite size of the 
survey volume. The SDSS data have been analyzed using both 
the p o wer spectrum method (see, e.g..lTegmark et alj|2004l; lllutsi 
2005 ; IPadmanabhan et all 120071 : iBlake et al.l 120071 : IPercival et al. 
2007, l20ld: lReidetal.1 l2of(if). and the correlation function 



Cabre & Gaztanag j 120081 ; iMartinez etail 120091 : ISanchez et ail 
2009 : lKazinetalj201Ch . 



method (see, e.g., lEisenstein et alj 1 20051 ; lOkumura et alj 1200 8t 
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The three major uncertainties while constructing a theoreti- 
cal prediction of the power spectrum or correlation function with a 
given cosmological model are the galaxy bias (the relationship be- 
tween galaxy and matter distributions), non-linear effects, and red- 
shift distortions. The knowledge of these uncertainties determines 
which analysis method and scale range we should use to obtain re- 
liable constraints on the dark energy and cosmological parameters. 

In this paper, we present the measurement of the spherically- 
averaged correla tion function from the SDSS DR7 lumino us red 
galaxy (LRG) teisenstein et al.l l200ll : lAbazaiian et alj|2009h sam- 
ple which provides a homogeneous galaxy sample and has the 
largest effective survey volum e to date for studying the linear 
regime teisenstein et alj l2005h . In Section [2] we introduce the 
galaxy sample and selection functions used in this study. In Sec- 
tion^ we describe the details of our method. In Sec. [4] we present 
our results. In Sec. [5] we check our results using systematic tests. 
We summarize and conclude in Sec. [6] 



© 0000 RAS 



2 Chuang, Wang, & Hemantha 



2 DATA 



The SDSS has observed one-quarter of the entire sky and per- 
formed a redshift survey of galaxies, quasars an d stars in five 
passbands u, q, r } i, and z w ith a 2.5m telescope dFukugita et al.l 
ll996tlGunn et al.ll998ll2006h . We use th e public catalog, the NYU 
Value-Added Galaxy Catalog (VAGC) jBlanton et al]l2005ri . de- 
rived fr om the SDSS II final public data release, Data Release 7 
(DR7) jAbazaiianet"ai]|2009h . We select our LRG sample from 
the NYU VAGC with the flag primTarget bit mask set as 32. K- 
corrections have been applied to the galaxies with a fiducial model 
(ACDM with fl m — 0.3 and h — 1), and the selected galaxies are 
required to ha ve rest-frame g-band abs olute magnitudes —23.2 < 
M g < -21.2 jBlanton & Roweisll2007l). The same selection crite- 
ria were used in previous work dZehavi et al. 2005 ; Eisenstei n et al, 
2005 1: lOkumura et alJl2008l ; ICabre & Gaztan aga 2008; Kazin et al 
2010). The sample we use is referred to as "DR7full" in l Kazin et al 
(2010). Our sample includes 87000 LRGs in the redshift range 
0.16-0.44. The average weighted redshift of the sample is 0.33. 

Spectra cannot be obtained for objects closer than 55 arcsec 
within a single spectroscopic tile due to the finite size of the fibers. 
To correct for these "collisions", the redshift of an object that failed 
to be measured would be a ssigned to be the s ame as the nearest 
successfully observed one lZeh avietai]|2005h . Both fiber colli- 
sion corrections and K-corrections have been made in NYU- VAGC 
jBlanton et alj2005l) . 

We have applied evolutionary correction based on the stellar 
synthesis model of an old and p assively evolving burst from z = 
10 by using the PEGASE code dFioc & Rocca- Volmerange|[T997h . 
The rest-frame g-band absolute magnitudes are passively e volved 
to z = 0.3 (see lZehavi et alj d2005l) : iEisenstein et all {2005)). 

We construct the radial selection function as a cubic spline fit 
to the observed number density histogram with the width Az = 
0.01 (see Fig. [j}. The NYU- VAGC provides the description of 
the geometry and completeness of the survey in terms of spheri- 
cal polygons. Although the c ompleteness of VAG C is determined 
based on the main galaxies dStrauss et al.1 12002). we adopt it as 
the angular selection function of our sample since the main galax- 
ies and LRGs should have similar angular selection functions (see 
the appendix of Zehavi et al. 2005). We drop the regions with 
comp leteness below 60% to avoid unobserved plates dZehavi et al.1 
2005). The Southern Galactic Cap region is also dropped. 



3 METHODOLOGY 

In this section, we describe the measurement of the correlation 
function from the observational data, construction of the theoretical 
prediction, and the likelihood analysis that leads to constraints on 
dark energy and cosmological parameters. We will also show that 
using one scaling parameter is sufficient for extracting information 
from the observed spherically averaged correlation function (see 
Sec. l3A3l 



3.1 Measuring the Two-point Correlation Function 

We calculate the comoving distances to every galaxy by assum- 
ing a fiducial model, ACDM with Q. m — 0.25. We use the two- 
point correlation funct ion estimator given by Landy and Szalay 
dLandv & Szalavll 19931) : 



DD(s) - 2DR(s) + RR(s) 
RRjs) 





Figure 1. The radial selection function of the LRG sample used in this 
study. The gray bars are computed from the sample and the blue line is 
the cubic spline fit of these bar values. We compute the radial selection 
function in the form of the number of galaxies per unit redshift instead of the 
number density in comoving coordinate, so that we don't need to assume a 
fiducial model while generating the random catalog with the radial selection 
function. 



where DD, DR, and RR represent the normalized data-data, data- 
random, and random-random pair counts respectively in a distance 
range. The bin size we use in this study is 5/i _1 Mpc. This estimator 
has minimal variance for a Poisson process. Random data should be 
generated according to the radial and angular selection functions 
of the data. One can reduce the shot noise due to random data by 
increasing the number of random data. The number of random data 
we use is 10 times that of the real data. While calculating the pair 
counts, we assign each data point a radial weight of 1/[1 + n(z) ■ 
P w ], where n(z ) is the radial selection function and P w = 4 ■ 10 4 
ft _3 Mpc 3 as in lEisenstein et al.1 d2005l) . The observed correlation 
function is shown in Fig. [2] 

3.2 Theoretical Two-Point Correlation Function 

We compute the linear power spectr a at z = 0.33 by using CAMB 
dLewis. Challinor. & Lasenbvl2000h . To include the damping effect 
of non-linear structure formation and peculiar velocities, we calcu- 
late the dewiggled power spectrum 



'dw{k) = P ;i „(fc)exp [ -7^2 



~\~Pnw (&) 



exp 



2kl 



where Pu n (k) is the linear power spectrum, P nw (k) is the no- 
wiggle or pure CDM p ower spectrum calculated with Eq.(29) in 
lEisenstein & Hul dl998l) . and is marginalized over with a flat 

prior over the range of 0.09 to 0. 120 

We then use the software package halofit dSmith et alJl2003h 
to compute the non-linear matter power spectrum: 



r halofit 



(AO 



halofit,i 



,(*) 



(fc) 



(3) 



(1) 



1 Although k+ can be computed by renorma lization perturbation theory 
i Crocce & Scoccimarro 2006; Matsubara 2007), doing so requires knowing 
the amplitude of the power spectrum, which is also marginalized over in 
this study. 
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Figure 2. The spherically-averaged two-point correlation function mea- 
sured from the SDSS DR7 data. The red triangles are the correlation func- 
tion computed wit h the LRG sa mple described in Sec. [2] The green circles 
are taken from lKazin et al.ll2010l in which the same fiducial model is used 
(ACDM with Q m = 0.25) but the bin size they use is lOfe^ Mpc. Our 
result shows excellent agreement with that of[Kazin et al. 2010. The black 
line is the average correlation function from LasDamas mock catalogs. The 
error bars are the square roots of the diagonal elements of the covariance 
matrix we have derived (see Sec. 13.31 . The violet dashed line is the mean 
model from our MCMC likelihood analysis (U m h 2 = 0.105, U),h 2 = 
0.0225, n s = 0.978, D v (0.35) = 1428Mpc). Note that an MCMC anal- 
ysis does not result in an accurate bestfit model <Lewis & Bridl e 2002). 



i{k) = Pd w {k)r ha i ofit {k), 



(4) 



where P^igfy nvj (k) is the power spectrum from applying halofit 
on the no-wiggle power spectrum and P„i(k) is the non-linear 
power spectrum. We compute the theoretical two-point correlation 
function by Fourier transforming the non-linear power spectrum. 
We show an example of the effect of applying dewiggle and halofit 
to the correlation function in Fig. [3] Clearly, the damping of the 
baryon acoustic oscillation (BAO) peak is accurately described by 
the dewiggled linear correlation function. Additional nonlinear ef- 
fects are only important on very small scales. 

The parameter set we use to compute the theoretical correla- 
tion function is {Dv(z), f2 m /i 2 , Q,bh 2 ,n s , fc*}, where Q m and fif, 
are the density fractions of matter and baryons, n s is the power law 
index of the primordial matter power spectrum, h is the dimension- 
less Hubble constant (Ho — lOO/i km s _1 Mpc _1 ), and Dv(z) is 
defined by 



Dv(z) 



(l + z) 2 D\ 



H(z) 



(5) 



where H(z) and Da(z) are the Hubble parameter and the angular 
diameter distance at the redshift, z. We set h = 0.7 while calculat- 
ing the non-linear power spectra. Assuming the true model of our 
universe is not far from a ACDM model, the dark energy and cur- 
vature dependence are absorbed by the effective distance, Dv(z). 
Thus we are able to extract constraints from data without assuming 
a dark energy model and cosmic curvature. 



Figure 3. An example of the effect of applying dewiggle and halofit to the 
correlation function. The black solid line is the linear correlation function 
without applying dewiggle and halofit yet. The red dotted line is the dewig- 
gled linear correlation function. The green dashed line is the dewiggled cor- 
relation function including nonlinear effects calculated using halofit. The 
damping of BAO is accurately described by the dewiggled linear correla- 
tion function. Additional nonlinear effects are only important on very small 
scales. 



3.3 Covariance Matrix 

We use the mock catalogs from the LasDamas simulation^ 
(McBride et al., in preparation) to estimate the covariance matrix 
of the observed correlation function. LasDamas provides mock cat- 
alogs matching SDSS main galaxy and LRG samples. We use the 
LRG mock catalogs from the LasDamas gamma release with the 
same cuts as the SDSS LRGfull sample, -23.2 < M g < -21.2 
and 0.16 < z < 0.44. We have diluted the mock catalogs to match 
the radial selection function of the observed data by randomly se- 
lecting the mock galaxies according to the number density of the 
data sample. We calculate the spherical-averaged correlation func- 
tions of the mock catalogs and construct the covariance matrix as 



N 



dm 



(6) 



where TV is the number of the mock catalogs, £ m is the mean of the 
m th bin of the mock correlation functions, and ^ is the value of 
m th bin of the k th mock correlation function. We test the normal- 
ity of the correlation functions from the LasDamas mock catalogs 
and find that they are well described by normal distributions (see 
Appendix IB!. 

The mock catalogs derived from N-body simulations require 
long computing times and are very limited in availability. It is inter- 
esting to investigate whether there is an easier, faster, and cheaper 
way to construct mock catalogs which could work as well as 
those derived from N-boday simulation. Tow ards this end, we have 
created 500 lognormal(LN) moc k catalogs jColes& Joneslll99ll: 
IPercival, Verde. & Peacock]|2004l) . and computed the spherically- 
averaged correlation functions from these. The details involved in 
creating LN mock catalogs are described in Appendix [A] We com- 
pare the correlation funcions from the LasDamas mock catalogs 



http://lss.phy.vanderbilt.edu/lasdamas/ 
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LasDamas mock catalogs 

LN mock catalogs 

input 2PCF of LN mock catalogs 





s(h Mpc) 



s (h Mpc) 



Figure 4. Spherical-averaged 2PCF of the mock catalogs. The black solid 
line is computed from the LasDamas mock catalogs. The red dashed line is 
computed from our lognormal mock catalogs. The error bars are the square 
roots of the dianogal elements of the covariance matrixes. The green dotted 
line is the input 2PCF for our lognormal mock catalogs. 



Figure 5. The normalized covariance matrix computed from 160 LasDamas 
mock catalogs. We show the covariance among 40 bins from < s < 
200/i _1 Mpc with bin size 5/i _1 Mpc. 



and LN mock catalogs in Fig. |4j the error bars indicate the square 
roots of the dianogal elements of the covariance matrixes. We also 
show the normalized covariance matrixes in Fig. [5] and [6] Clearly, 
the results from the LasDamas mocks and our LN mocks are very 
similar to each other. In particular, the input correlation function is 
accurately recovered by analyzing the LN mock catalogs. Note that 
the LN mocks give larger errors on all scales, and on scales smaller 
than ~ 40/2 _1 Mpc, the LN mock catalogs give much larger errors 
than the LasDamas mock catalogs (see Fig.[7j>. 

We use the covariance matrix computed from the LasDamas 
SDSS mock catalogs, since these are more realistic than the log- 
normal mock catalogs, and give smaller errors for the measured 
correlation function. It is interesting to note that in the absence of 
mock catalogs derived from cosmological N-body simulations, log- 
normal catalogs can be used for a conservative estimate of the co- 
variance matrix of the correlation function. 



1 . 



-0.5 





Figure 6. The normalized covariance matrix computed from 500 lognormal 
mock catalogs. We show the covariance among 40 bins between the scale 
range, < s < 200/t~ 1 Mpc, with the bin size, 5h. _1 Mpc. 



3.4 Likelihood 

The likelihood is taken to be proportional to exp(— x 2 /2), and \ 2 
is given by 

X = ^2 [£th(Si) ~ £o6s(Si)] Cy 1 [ZthiSj) - ^obs(Sj)] (7) 
i,i=l 

where Nuns is the number of bins used, ^ t h is the theo- 
retical correlation function of a model, and £ bs is the ob- 
served correlation function. Note that (,th(si) depends on 
{Dv-(z), Qmh 2 , £l b h 2 , n 3 , fc*}. 

In principle, we should recalculate the observed correlation 
function while computing the \ 2 f° r different models. However, 
since we don't consider the entire scale range of the correlation 
function (we only consider s = 40 — 120 /i~ 1 Mpc in this study), 
we might include or exclude different data pairs for different mod- 
els which would render \ 2 values arbitrary. Therefore, instead of 
recalculating the observed correlation function, we apply the in- 
verse operation to the theoretical correlation function to move the 



parameter dependence from the data to the model, thus preserving 
the number of galaxy pairs used in the likelihood analysis. 

Let us define T as the operator converting the measured cor- 
relation function from the fiducial model to another model, i.e., 

Us{s) = miitis)). (8) 

where ^fL(s) is the observed correlation function assuming the 
fiducial model. This allows us to rewrite \ 2 as 

X 2 EE { T " 1 [^(s i )]-$ / "(^)}C7i ir 

' {T~ X \£th{sj)] — fobs(sj)} j (9) 

where we have used Eqs.(6jl and (§}. 

To find the operator T, note that the fiducial model is only 
used in converting redshifts into distances for the galaxies in our 
data sample. In the analysis of galaxy clustering, we only need the 
separation of a galaxy pair, and not the absolute distances to the 
galaxies. For a thin redshift shell, we can convert the separation of 
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LN mock catalogs 
LasDamas mock caialojis 



s(h Mpc) 



s (h Mpc) 



Figure 7. Square roots of the dianogal elements of the covariance matrixes. 
The black dashed line is computed from the LasDamas mock catalogs. The 
red solid line is computed from lognormal mock catalogs we create. One 
can see that LN mock catalogs have larger covariance at smaller scale and 
two line are close at larger scale. 



one pair of galaxies from the fid ucial model to another m odel by 
performing the scaling (see, e.g. JSeo & Eise nstein (200^)) 



Hf id (z 

i — 

H(z) 



D A (z) 



D 



fid/ 



(10) 



where 9 is the angle between the radial direction and the direction 
of th e line connecting the pa ir of galaxies. 

lEisenstein et al. ( 2005) argued that we can use one rescal- 
ing parameter, Dy{z), to convert the observed correlation function 
from the fiducial model to another model as long as the new model 
is not very different from the fiducial one, and the redshift range of 
the sample is not large. Then the separation of one pair of galaxies 
is converted from the fiducial model to another by 



Dyjz) 
D f v * d (z) 



(11) 



In this section, we discuss methods with one and two rescaling 
parameters, and show that these two methods are equivalent for 
spherically-averaged data when certain conditions hold (see Sec. 
[3A3J. 



3.4.1 Using One Rescaling Parameter 

From eq. ( Hit , the observed correlation function with the different 
model can be written as follows: 



frd ( Dy d (Zeff) 



Dv(Zeff) 



(12) 



where z e ff is the effective redshift of the sample and Dyiz) is 
defined by Eq.(|5}. 

The effective redshift we use in this study is z e ff = 0.33. 
Since the results are insensitive to z B f / (see Sec. [5), we rescale our 
result to z e f / = 0.35 for comparing with previous works. Eq. d!2t 
can be rewritten as 



\D J V (Zeff) ) 



(13) 



We can apply the same inverse rescaling operation to the the- 
oretical correlation function: 



T (ftft(s)) - &h —j-j- 

\ D V ( Z eff) 



X 2 can be calculated by substituting eq. d 14fc into eq. J5J. 



(14) 



3.4.2 Using Two Rescaling Parameters 

From eq. <10t , we can convert the spherically-averaged correlation 
function from some model to the fiducial model by 



fid/ 



d9w(s,9)x 




H{z) 



+ ssinl 



, Da(z) 
D{ td (z) 



where the weighting function w(r, 6) is given by 

n DD (s, 6) 



J* d6n DD (s, 9) ' 



(15) 



(16) 



where udd(s, 9) is the number density of the data pairs. We define 
inverse operation, T~\ directly since T is not necessary in our 
calculation. We now apply the inverse operation to the theoretical 
correlation function: 

r _1 (6h(s))= / d9w(s,9)x 



Hf id (z) 



Da(z) 
D% d (z) 



(17) 



X 2 can be calculated by substituting eq. J17t into eq. {9]l. 



3.4.3 Equivalence of Using One and Two Rescaling Parameters 
for Spherically-Averaged Data 

We now show that using one and two rescaling parameters while 
calculating the spherically-averaged correlation function are equiv- 
alent to first order in approximation. If the size of the survey is 
much larger than the scales of interest, udd(s,6) would be pro- 
portional to s sin 9. Hence 



w(s,6) 



s sm f 



/»7T 

Jo 



s sin9d9 



sin 9 
2 



(18) 



Next, if the model is close to the fiducial model, we can just con- 
sider the first order terms of D v /D{; d , H fld /H, and D A /D f A id 
which can be written as following: 



Dy 
nf id 



1 + 5 V 



H fic 



1 + Sr 



Da 

A 



1 + S a , 



(19) 



where \8v \ , \S r | , \8 a | 1. From the definition of Dv (see Eq.l5l). 
one can obtain a simple relation, 35v — S r + 2S a . Let's consider a 
power law correlation function: 



(20) 
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where p is a real number. Eg.d 17t can be rewritten as 



de—- I J [s cos 0(1 + 5 r )T + [s sin 0(1 + 5 a )]' 



dd sin [cos 2 0(l + 25 r ) + sin 2 0(1 + 28 a ) 



o 

1 + p 



dO sin Tl + p(cos 2 65 r + sin 2 88 a 
8 r + 28 a 



s p (l+p6 v ) 
Dv 

pifid ' 



(21) 



The proof can be generalized to any function which can be ex- 
pressed as 



U(s) 



+ . 



(22) 



where pi are real numbers. 

To measure the spherically-averaged correlation function, we 
have shown that using one rescaling parameter, Dv, and two 
rescaling parameters, H and Da, are equivalent as long as the 
scales of interest are relatively small compared to the survey length 
scale, and the constraint on Dv is tight enough. A similar statement 
can be made for the spherically-averaged power spectrum analysis. 
We have verified that these two rescaling methods give similar re- 
sults. 



3.5 Markov Chain Monte-Carlo Likelihood Analysis 

We use CosmoMC dLewis & Bridle! |2002|) in a Markov Chain 
Monte-Carlo likelihood analysis. The main parameter space that 
we explore is {Q m h 2 , ftbh 2 , n B , Dv(z e ff), fe*} and the prior 
ranges are {(0.025,0.3), (0.01859,0.02657), (0.865,1.059), 
(725,1345), (0.09,0.13)} respectively. The dependence on h, 
the curvature, and dark energy parameters are absorbed into 
D v (z ef f). 

We marginalize over the amplitude of the correlation func- 
tion; this is equivalent to marginalizing over galaxy biasxag x 
where ag is the matter power spectrum normalization parameter 
and rp is the linear ratio between the correlation function in the 
redshift space and real space which can be derived from the linear 
redshift distortion parameter dKaiseil[l987l) . Since the LRG data 
alone cannot give tight constraints on Q, n h 2 and n s , we apply flat 
priors (+7awMAp) on them which are wide enough so that CMB 
constraints will not be double counted. In other words, the effect 
from the wide flat priors could be ignored when combining our fi- 
nal results with CMB data. We also marginalized over k, over the 
range of 0.09 to 0.13 (see Sec.3.2). 



4 RESULTS 

In this section, we present the model independent measurements of 
the parameters we explore, {Dv(0.35), Q. m h 2 }, obtained by using 
the method described in previous sections. Although, the effective 
redshift we use is 0.33, the average weighted redshift, we rescale 



all our results to z e f / 
easily by 

£>v(0.35) = IM0.33) 



0.35 for comparing with previous work 



D VJid (0M) 



1.054D V (0.33). (23) 



£>v,/ ld (0.33) 

We have checked that the results is insensitive to the effective red- 
shift in Sec. [5] 

We derive the model independent measurements of H and Da 
for comparison with 2D results. We also apply the method on two 
subsamples (two redshift slices) as a systematic test. 

We validate our method by applying it to the LasDamas mock 
catalogs, and find that our measurements are consistent with the 
input parameters of the simulations. 

We derive constraints on dark energy and cosmological pa- 
rameters by combining our res ults with other d ata sets including 
WMA P7 dKomatsu etalB oiO) and Union2 SN iAmanullah et alj 
l2010h . 

Finally, we compare our results with previous works. 

4.1 Model Independent Constraints on Dv (0.35) 

Without assuming a dark energy model or a flat Universe, we 
find that D v (0.35) = 14281^ Mpc and r s (z d )/D v (0.35) = 
0.1143 ± 0.0030, where r a (z a ) is the comovin g sound horizon 
at the drag epoch calculated with the eq. (6) in Eis enstein & Hul 
dl998l) . 

Fig. [8] shows one and two-dimensional marginalized con- 
tours of the parameters, {/JV(0.35), Q, m h 2 , r B (zd)/ Dv(0-35), 
A(0.35)}, where 



A(z) 



D 



v(z 



(24) 



The measurements and the covariance matrix are listed in Table 
[T] and [2] The bestfit model from the MCMC likelihood analysis 
has \ 2 ~ 6.32 for 16 bins of data used (in the scale range of 
40/i _1 Mpc< s < 120ft _1 Mpc with the bin size = 5/i _1 Mpc), 
for a set of 6 parameters (including the overall amplitude of the 
correlation function). 

The scale range of the correlation function we have selected 
is s — 40 — 120 /i~ 1 Mpc. In this range, the scale dependence 
of the redshift distortion and galaxy bias is small. We cut the tail 
of the correlation function at s — 120 /i _1 Mpc because the high 
tail (large correlation at large scales) cannot be fitted to any con- 
ventional model, and could be due to systematic error or sample 
variance (see further discussion in Sec.|5j. 

At this point, we assume the high tail is simply due to sample 
variance, and might disappear when much larger data sets become 
available. Unlike previous analyses by other groups, we apply very 
weak flat priors (+7ctwmap7) on Qth 2 and n a instead of fixing 
them to the bestfit values from CMB data. 



4.2 Model independent measurements of //(0.35) and 
Da(0-35) 

In this section, we apply the method with two scal- 
ing parameters described in sec 13.4.21 to measure H 
and Da- The main parameter space that we explore is 
{Q m h 2 ,Q b h 2 ,n s ,H(z e ff),DA(z e ff),k i ,} and the prior ranges 
are {(0.025,0.3), (0.01859,0.02657), (0.865,1.059), (41,123), 
(723,1343), (0.09,0.13)} respectively. We obtain the model 
independent measurements, i/(0.35) = 83^^ km s _1 Mpc _1 and 
L>a (0.35) = 10891^ Mpc, from the LRG data alone (see Table 
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20014001600 



a 0.1 

0.08 




120014001600 



0.08 0.12 0.16 





mean 


a 


lower 


upper 


#(0.35) (km s^Mpc" 1 ) 


83 


17 


68 


96 


D a (0 35) (Mnc) 


1089 


93 


1002 


1182 




0.105 


0.017 


0.089 


0.122 


#(0.35) *r s (z d ) 


13500 


2700 


11200 


15500 


r s (z d )/D A (0.35) 


0.151 


0.012 


0.140 


0.161 


A(0.35) 


0.432 


0.026 


0.408 


0.457 



120014001600 



0.08 0.12 0.16 




0.110.120.13 




Figure 8. 2D marginalized contours (68% and 95% C.L.) for D v (0.35), 
Qmh 2 , r s (z d )/ Dy (0.35), and A(0.35). The diagonal panels represent 
the marginalized probabilities. 



Table 3. Measured cosmological parameters with flat prior 0.01859 < 
fl b h 2 < 0.02657, 0.865 < n s < 1.059 (±7a WMAP7 ), and 0.09 < 
fc„ < 0.13/i _1 /Mpc. The standard deviations and the margenalized 
bounds (68%) are listed as well. There are three derived measurements, 
#(0.35) * r 3 (z d ), r s (z d )/D A (Q.35), and A(0.35), in the table. 





mean 




lower 


upper 


ZV(0.35)(Mpc) 


1428 


74 


1355 


1502 


n m h 2 


0.105 


0.016 


0.090 


0.121 


r s (z d )/D v (0.35) 


0.1143 


0.0033 


0.1113 


0.1173 


A(0.35) 


0.439 


0.020 


0.419 


0.459 



Table 1. Measured cosmological parameters with flat prior 0.01859 < 
fl b h 2 < 0.02657, 0.865 < n a < 1.059 (±7crwMAPr), and 0.09 < 
k t < 0.13/i — 1 /Mpc. The standard deviations and the marginalized 
bounds (68%) are listed as well. There are two derived measurements, 
r a (z d ) /D v (0.35) and A(0.35), in the table. 




A 



Figure 9. 2D marginalized contours for 68% and 95% for H(z = 0.35), 
D A (z = 0.35), Q m h 2 , #(0.35) * r s (z d ), r s (z d )/D A (0.35), and 
A(0.35). The diagonal panels represent the marginalized probabilities. 



0. Table[4]shows the normalized covariance matrix of {#(0.35), 
Da (0.35), n m h 2 , H(0.35)*r 3 (z d ), r s (z d )/D A (0.35),A(0.35)}, 
and Fig. |9] shows the 2D marginalized contours of this parameter 
set. 

Although using two rescaling parameters on the spherically- 
averaged correlation function cannot give better constraints on the 
cosmological parameters, it gives the model independent measure- 
ments of H and D A which cannot be derived directly from the 
measurement of Dv ■ These can be compared to our result for the 
two-d imensional two-point correlation function ( Chuang & Wang 
l201ll) . #(0.35) = 82.1±*;9 kms -1 Mpc -1 and £> A (0.35) = 



£V(0.35) tt m h 2 



D v (0.35) 



A(0.35) 



£V(0.35) 


1 


-0.7890 


-0.5561 


-0.1727 




-0.7890 


1 


0.0056 


0.7305 


£V(0.35) 

A(0.35) 


-0.5561 


00056 


1 


-0.6181 


-0.1727 


0.7305 


-0.6181 


1 



Table 2. Normalized covariance matrix with flat prior 0.01859 <Q. b h 2 < 
0.02657, 0.865 < n s < 1.059 (±7ct WM ap7), and 0.09 < fc* < 

o.m-VMpc- 



1048lgg Mpc. Not surprisingly, information is lost in the spheri- 
cal averaging of data. 



4.3 Model independent measurements of r B (z d ) /Dv {z) from 
two redshift slices 

To explore the redshfit dependency of the measurements. We apply 
the method of one rescaling parameter on two subsamples have z = 
0.16 — 0.36 and z — 0.28 — 0.44. The average weighted redshifts of 
these two samples are 0.28 and 0.36 respectively. Both subsamples 
have ~ 60000 galaxies. We find r s (z d )/D v (0.28) = 0.141 ± 
0.012 and r s (z d )/D v (0.36) = 0.1146 ± 0.0068. Because of the 
overlapping of the samples, there is covariance between these two 
measurements. To estimate the covariance, we apply the method on 
40 mock catalogs (indexed from Ola to 40a) and find the correlation 
coefficient r = 0.20. To combine these measurements with other 
data sets, one should add the following term to the \ 2 '■ 



XLRG2Z = &.LRG2Z 



6688 -2437 
-2437 22200 



where 



r s {z d )/D A (0.28) -0.1413 
r s {z d )/D A (0.36) -0.1146 



(25) 



(26) 
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#(0.35) 


£> A (0.35) 




#(0.35) *r s (z d ) 


r s (z d )/D A (0.35) 


A(0.35) 


H(0.35) 


1 


0.4028 


0.2133 


0.9744 


-0.5733 


-0.3460 


D j (0.35) 


0.4028 


1 


-0.4829 


0.5300 


-0.8440 


-0.1210 




0.2133 


-0.4829 


1 


0.0025 


-0.0139 


0.6035 


H(0.35)*r s {z d ) 


0.9744 


0.5300 


0.0025 


1 


-0.5861 


-0.4723 


r s (z d )/D A (0.35) 


-0.5733 


-0.8440 


-0.0139 


-0.5861 


1 


-0.2033 


A(Q.35) 


-0.3460 


0.0908 


0.6035 


-0.4723 


-0.2033 


1 



Table 4. Normalized covariance matrix with flat prior 0.01859 < Q, b h 2 < 0.02657, 0.865 < n a < 1.059 (±7a WM APT), and 0.09 < fc» < 

o.m^VMpc- 



The constraints of cosmological parameters from combining 
CMB and SNe data set for the owCDM model are shown in Ta- 
ble. [5] The measurements are similar to using the whole sample 
(z=0. 16-0.44) but the constraints are worse. It might be due to some 
tension between the measurements of these two subsample which 
might be caused by higher noise level from smaller samples or by 
the evolution of the dark energy which would need a more com- 
plex model than owCDM. We will explore this issue in our future 
research. 



4.4 Constraints on owCDM Model 

We now present the cosmological parameter constraints for the 
owCDM model (non-flat Universe with a constant dark energy 
equation of state). Table. [5]also shows the constraints from cosmo- 
logical microwave background (WMAP7) and supernova (Union2 
compilation) data and their combination with SPSS LRG dat a. To 
include the constraints from WMAP7 dKomatsu etail fcoiO). wc 
us e the constraints on the CMB shift p arameters {R,l a } and z* 
bv lWang. Chuang. & Mukheried d201lh (see Appendix [C}. To cal- 
culate the constraints from Union2 SNe, we use the add-on code 
for cosmoMC which can be download from the website of Union2 
SNcB For a given model, one could obtain \ 2 f° r each data set, i.e. 
Xcmb an d Xsivfl To include the constraint we obtained from the 
galaxy clustering data, one should add the following term to the \ 2 
with 



Xlrg — 



r a (z d )/D v (0.35) -0.1143 



0.0033 



(27) 



Combining all three data sets, LRG, CMB, and SNe, and as- 
suming the owCDM model, we find that Sl k = -0.0032±g;^, 
and w = -1.010 + n ?,^, which is consistent with ACDM model 



(in agreement with previous work, see e.g..lSerra et alj 2009; Wang 
l2009tlMortonson. Hu. & HuteretfeOldlzhao & Zhang||201ol) . FieT 
[To] compares the constraints on w and Qk in the owCDM model. 
We can see that the addition of SDSS LRG data significantly tight- 
ens the constraints on dark energy and cosmological parameters. 



4.5 Validation Using Mock Catalogs 

In order to validate our method, we have applied it to the 2D 2PCF 
of 40 LasDamas mock catalogs (which are indexed with 01a-40a). 
Again, we apply the flat and wide priors (±7<twmap7) on fl b h 2 



3 http://supernova.lbl.gov/Union/ 

4 While computing X 2 SN , we use the covariance matrix with systematics 
to obtain more reliable constraints from SNe 




-0.5 



Figure 10. 2D marginalized contours for 68% and 95% for w and 
(owCDM model assumed) from WMAP7 (dashed blue), WMAP7+Union2 
SN (dotted red), and WMAP7+Union2 SN+LRG (solid black). The straight 
solid black lines indicate that w = — 1 and = 0. 





mean 


<7 


input value 


£V(0.35)(Mpc) 


1349 


69 


1356 


Q m h 


0.120 


0.015 


0.1225 


r s (z d )/D v (0.35) 


0.1205 


0.0059 


0.1175 


A(0.35) 


0.441 


0.026 


0.452 



Table 6. The mean, standard deviation, and the 68% C.L. bounds of 
{IV(0.26), Qmh 2 , r s (2 d )/ZV(0.35), A(0.35)} from the LRGfull 
mock catalogs of the LasDamas simulations. Our measurements are con- 
sistent with the input values within ltr. 



and n s , centered on the input values of the simulation (Qbh 2 — 
0.0196 and n s = 1). 

Table [6] shows our measurements of {Z?v(0.35), Q. m h? , 
r 3 (z d ) / Dy (0.35), ^4(0. 35)} from the LasDamas mock catalogs 
of the SDSS LRG sample. These are consistent with the input pa- 
rameters, establishing the validity of our method. 



4.6 Comparison with other analyses 

There have been several analyses of the cluste r ing o f the 
SDSS LRG spectroscopic sample. IPercival et alj d201 Ol) use 
power spectrum analysis to extract the BAO signals. The sam- 
ple they use include not only LRGs but also main galax- 
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owCDM 





f2 m 


fix 


w 


Ho 


n m h 2 




CMB 
CMB+SN 
CMB+SN+LRG 
CMB+SN+LRG2Z 


n „, ,-+0.072 
u ""°-0.077 
0.348i«- 

0.278 ± 0.032 


67+ 073 
-0.071 

676+ 052 
-0.04S 

0.725 ± 0.018 
0.725 ± 0.026 


-1 04+ " 62 
" -0.63 

-1.11 ± 0.11 

-1 Ol+O-O 4 ^ 
-0.045 

-1 008+ 053 
-0.0S4 


65. 7 ;?- 2 2 

63.4+|"| 
69.9 ± 2.3 
70.3 ± 4.2 


13S6+ 0059 
U.13&b_ 006Q 

1361+ " 0061 
O.13bl_ 0060 

n -iqe 7 +0.0061 
136 + 0061 

"'"''-o.ooea 


006+ " 058 
""""-0.048 

-0 025+ " 022 
U - U2;> _0.019 

-0 0032+ 0074 
U.OU3-2_ 0072 

-0.003 ± 0.011 



Table 5. Constraints of the cosmological parameters from various data combinations with owCDM model assumed, where LRG is using the fiducial result of 
this paper (eq. [27} and LRG2Z is using the measurements from two subsamples as described in Sec. 14.31 There are two inferred parameters, f2 m and Q^, in 
this table. 



ies sample. They obtain the c onstra int on r a /Dv (0.35) = 

0. 1097 ± 0.0036. iReid et al.l d2010h apply power spectrum 
analysis on the reconstructed halo density field derived from 
SDSS DR7full (flux-limited LRG) sample. Their measurement 

of r s /Dy (0M) = .1097±g;go42 = °- 2 ) which is 

similar to iPercival et al. I d2010h . There is about lcr difference 
between their results and ours. However, in IReid et al.l feOlOl) . 
they show a dependence on the minimum scales of the range, 

1. e. r s /LV(0.35) = 0.11181™^ (fe mo3: = 0.15) and 
r s /D v (0-35) = 0.1136t° :l°r° 0™^ = 0.1). Their results are 
more consistent with ours when a smaller k max (i.e., a larger min- 
imum scale) is used, which represents a more conservative choice 
for the scale range. 



constant. In other words, they assumed a constant shift from some 
systematic error. Although this systematic error is unknown, we 
could minimize its effect by using smaller scale. The reason is that 
the correlation function has larger value at smaller scale so that the 
results are less sensitive to the shift. We choose s = 120h~ Mpc 
as our boundary for the large scale and show that the results are 
insensitive to the constant shift by lowering down the data points 
of the observed correlation function by 0.002. We find that fl m h 2 
varies by la and r 3 /Dy (0.35) only varies by 0.2a. Therefore, our 
measurement of r s /Dv (0.35) is robust to the systematic shift. This 
is another indication that our measurement of r s (zd) / Dy (0.35) 
(but not that of Q. m h ) is robust. 



5 SYSTEMATIC TESTS 

Table. [7] shows the systematic tests that we have done varying key 
assumptions made in our analysis. These include the range of scales 
used to calculate the correlation function, the nonlinear damping 
scale, an overall shift in the measured correlation function due to a 
systematic error. 

We vary the effective redshift (from z e ff = 0.33 to 2 e // = 
0.35) used to calculate the theoretical model. We rescale the results 
to z — 0.35 for comparison and find that the results are insensitive 
to the effective redshift. 

We also test the sensitivity of our results to the nonlinear 
damping scale, fc*. Although k+ can be predicted accurately in 
the real space iCrocce & Scoccimarrd[200q ; iMat subara 20071), in 
the redshift space, it would also depend on the redshift distortions 
which cannot be well determined from the spherically-averaged 
correlation function. In table. [7j one can tell that the results are 
not sensitive to fc*. 

In principle, the range of scales chosen for the analysis should 
be as large as possible, in order to derive the tightest constraints. 
However, we do not use the small scales (s < 40 h~ Mpc), where 
the scale dependence of redshift distortion and galaxy bias are not 
negligible and c annot be accura t ely de termined at present. Accord- 
ing to Fig. 5 in Eisenstein et al. (2005), these effects are negligible 
at s > 40 fa -1 Mpc. We vary the minimum scales used and find 
that the r s (zj) / Dv (0.35) is insensitive to it but VL m h 2 is not. This 
indicates the robustness of the measurement of r a (zd) / Dv (0.35) 
(but not Q, m h 2 ) from this paper. 

On larger scales (s > 130 ft _1 Mpc), the observed correla- 
tion is significantly higher than expected in conventional models 
of galaxy clustering. This high tail problem was reported i n pre- 
vious work, see, e.g., Eisenstein et al. (2005). lHutsil J2005T) . and 
I Sanchez eTa l. (2009). They found that the observed correlation 
function could be fitted better by lowering all the data points by a 



6 CONCLUSION 

We have presented our first results for the model independent con- 
straints on dark energy from the spherically-averaged correlation 
function of SDSS DR7 data, using an MCMC likelihood anal- 
ysis. Our constraints on {£V(0.35), fi m /i 2 , r s (zd)/ Dv(0-35), 
A(0.35)} are summarized by Table Q] and [2] Applying these re- 
sults to constraining a constant dark energy equation of state with- 
out assuming a flat Universe (the owCDM model), and combin- 
ing with WMAP7 and Union2 SN data sets, we find that fi*. = 
-0.0032±g'°° 7 4 and w = -l.OlOtg;^, consistent with a flat 
universe with a cosmological constant (Qk = 0, w = — 1). 

We have also measured the model independent constraints 
of H(0.35) and Da(0.35) from the spherically-averaged correla- 
tion function from SDSS DR7 LRGs, as a baseline for compar- 
ison with constraints from studies of the 2D correlation function 
(see lChuang & Wand fcOl lh ). We find that {#(0.35), D A (0-35)} 
from the spherically-averaged correlation function provide much 
weaker constraints than the 2D correlation function; this is as ex- 
pected since spherically-averaging reduces the amount of informa- 
tion extractable from data. 

The correlation function analysis is expected to be a more 
robust way to extract the BAO signals than the power spec- 
trum analysis, because one can easily get rid of the system- 
atic uncertainties such as the redshift distortion, the galaxy bias, 
and the non-linear effect by cutti ng off the small scale range 
dSanchez]~B augh. and Anguloll2Q08l) . 

The power of the correlation function analysis is limited at 
present by the available data. The correlation funciton that we have 
measured from the SDSS DR7 data has a high tail (larger than ex- 
pected correlations) at large scales (s > 120) (see Fig.[2j. Whether 
this high tail is simply due to the sample variance or some other 
systematic issue, e.g., extinction correction, will only become clear 
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D v (0.35)(Mpc) 



Q m h 2 r s (z d )/D v (0.35) 



A(0.35) 



fiducial model 


1428 


z.ft = 0.35 


1427 


fc* = 0.11 


1426; 


s=20-120 


1398; 


s=60-120 


1418 


s=40-100 


1448; 


s=40-140 


1393 


shift = 0.002 


1388 



+ T4 
-73 
+72 
-71 
+73 
-72 
+64 
-62 
+94 
-93 
+96 
-89 
+90 
-87 
+80 
-77 



0.105 ±0.016 
105+ - 016 
0.106 ±0.015 
0.116 ±0.012 
0.108 ±0.025 

n 114+°018 
ull ^-0.019 

0.110 ±0.018 
0.122 ±0.020 



1143+00029 
u.±±to_ 0031 

H44+0.0030 
u.±-L<i<i_ Q031 

n 1144+ 0030 
u.ii44_ 0031 

0.1139 ±0.0028 

1149+0 0030 
u.i±±ti_ 0Q33 

n 1,1+0.0050 
u - lxl -0 .0058 

1164+ 0040 

u ' ilu *-0 .0047 
H36+ 0032 



0.439 ± 0.02 
0.439 ±0.02 
0.439 ±0.02 
0.453 ±0.015 
0.438 ±0.03 

n 463+0-033 
u,wo -0 .032 
4 o 7 +0.025 
U -^ J ' -0.023 

0.459 ± 0.023 



Table 7. This table shows the systematic tests with the scale range, the fiducial model used, the effective redshift, the damping factor, and the shift from a 
systematic error. The fiducial results is obtained by assuming ACDM model with f2 m = 0.25 as fiducial model considering the scale range (s = 40 — 
120 /i _1 Mpc), using the effective redshift (z e f j = 0.33)), and the damping factor, k t , marginalized over with the a flat prior (0.09 < fc* < 0.13 /iMpc~ 1 ). 
The other results are calculated with only one quantity different from the fiducial one. n s = 0.963 and Q b h 2 = 0.02258 are marginalized with the same flat 
priors (±7ct\ymap7) m this paper. 



as more ambitious galaxy survey data become available in the fu- 
ture (e.g., from BOSSlI or Euclid^). 
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APPENDIX A: LN MOCK CATALOGS 

One convenient way to generate mock galaxy catalogs for cal- 
culating covariance matrix is using lognormal random fields 
which can approxima te the present-day non-linear fluctuation field 



Coles & Jones 1991). We create 500 lognormal (LN) de nsity fields 
Coles & JoneJl99ll ; |Percival, Verde, & Peacock 2004) on a 512 3 



grid with box length 4096 /i _1 Mpc. We then draw a random Pois- 
son variable with mean given by the selection functions and lognor- 
m al field to create the mock ca t alogs. We follow the steps described 
in lPercival. Verde. & Peacockl d2004l) except that we don't cut the 
input power at 0.25 Nyquist frequency because it makes the re- 
stored correlation function deviate from the input one. With a input 
correlation function, £(r), the Gaussian field correlation function is 
obtained by 



te(r)=ln[l+£(r)], 



(Al) 



and this can be Fourier transformed to the power spectrum, Pa(k). 
A Gaussian density field So (r) is generated on the grid with this 
power spectrum, and the corresponding lognormal field is calcu- 
lated by 




Normal order statistic means 



Figure Bl. The normality probability plot of the first bins, 40 < s < 
45h~ 1 Mpc, of correlation functions from the LasDamas mock catalogs 
we use for estimating the covariance matrix. The vertical axis is from the 
ordered values of the first bins from 160 mock catalogs which have been 
shifted and linear rescaled to have zero mean and unity variance. The hor- 
izontal axis is from the expected values of 160 ordered gaussian random 
numbers with zero mean and unity variance. That it is an approximate 
straight line means the bin values are approximately normally distributed. 



SLN(r) 



exp 



Mr) - f 



(A2) 



where 1 + Sln (r) is the lognormal density field which is always 
positive by definition and <Tq is the variance of the Gaussian density 
field which can be calculated by 



2 

0"G 



Po 



(k x . + k y + k z ) ; 



(A3) 



where N, 



2tt 
L 



grid 
N n „,„ 



s the number of grid points, k mn = 

, L is the box length, and m — x, y, or z. 

Then, the mock catalogs can be constructed by drawing the 
Poisson random variables with the means given by this lognormal 
field and the selection function of the galaxy survey. 

To compute the correlation function of these mock catalogs, 
one should create the random data on the same grid as well to can- 
cel out the effect of the finite size of the grid. The input correlation 
function in this study is the theoretical correlation function with 
parameters (f2 m = 0.25, fit = 0.04, h — 0.7, n s — 1) which 
are the same as the input parameters of the LasDamas simulations. 
We fix k* = 0.11 and the amplitude is adjusted to fit the averaged 
correlation function from the LasDamas mock catalogs we use. We 
are not fiting the observed correlation function because we want to 
find out is whether the LN mock catalogs could behave as good as 
LasDamas mock catalogs while estimating the covariance matrix. 



APPENDIX B: NORMALITY TEST 

We check the normality of the correlation functions from the Las- 
Damas mock catalogs by showing the normal probability plots of 
the first and last bin we use (see Fig. lBll and lB2t . One can tell that 
they are well described by a normal distribution. 
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Normal order statistic means 



Figure B2. The normality probability plot of the first bins, 115 < s < 
120/i -1 Mpc, of correlation functions from the LasDamas mock catalogs 
we use for estimating the covariance matrix. The vertical axis is from the 
ordered values of the first bins from 160 mock catalogs which have been 
shifted and linear rescaled to have zero mean and unity variance. The hor- 
izontal axis is from the expected values of 160 ordered gaussian random 
numbers with zero mean unity variance. That it is an approximate straight 
line means the bin values are approximately normally distributed. 



APPENDIX C: CMB DISTANCE PRIORS 

Wan g & Mukheried J2007t) showed that CMB shift parameters 
{l a , R), together with fl b h 2 , provide an efficient and intuitive sum- 
mary of CMB data as far as dark energy constraints are concerned. 
It is equivalent to repl ace Q b h 2 with z f , the redshift to the photon- 
decoupling surface (Wang 20091) . 

The CMB shift parameters are defined as (Wang & Mukherjed 

120071) : 



R = y/Q m Hl r(z»), l a = irr(z*)/r s (z*), 



(CI) 



and Zi, is the re dshift to the photon-deco upling surface given by the 
fitting formula fHu and Sugivama|[l996l) : 

z, = 1048 [l + 0.00124(fi(,/T, 2 )" ' 738 ] [l + 9l (n m h 2 ) 92 ] , (C2) 
where 

,2\-0.238 



9i = 



92 = 



0.0783 (Q b h 2 



l + 39.5(fi„/i 2 ) 0763 
0.560 



l + 21.1(fi 6 /i 2 ) 1 - 81 
The comoving distance to an object at redshift z is given by: 

r(z) = cH.7 1 |^r 1/2 sinn[jfi fc | 1/2 r( Z )], 



r(*) = 



|"k| 
z dz' 

E ( z ') 



(C3) 
(C4) 

(C5) 



E(z) = H(z)/H 



where sinn(x) = sin(x), x, sinh(a;) for fi^ < 0, fi& = 0, and 
fife > respectively; and the expansion rate the universe H{z) is 
given by 

# 2 M=(£) 2 (C6) 

= H% [fi m (l + zf + fi r (l + zf + fi fc (l + z) 2 + QxX(z)] , 

where fi m + fi r + fife + fix = 1, and the dark energy density 
function X(z) is defined as 



X(z) 



px(z) 

px(oy 



(C7) 



Note that fi r <C fim, thus the fi r term is usually omitted in dark 
energy studies, since dark energy should only be important at late 
times. 

The comoving sound horizon at redshift z is given by 

H c s dt' 



r s (z) = 



dz' ■ 



E(z'Y 



cm 



da' 



\/3(l + R b a')a' i E 2 {z') 



(C8) 



where a is the cosmic scale factor, a = 1/(1 + 2), and a 4 E 2 (z) — 
fi m (a + a oq ) + fife-a 2 + Q. x X{z)a i , with a cq = fi rad /fim 
I ;1 + z cq ), and 2 cq = 2 5 x 10 4 fi m fe 2 (Tb M a/2.7 K)' 4 . The 
sound speed is c s = l/y/ 3(1 + Rb a), with R b a = 3pb/(4p 7 ), 
R~ b = 31500fi D /i 2 (r C A/B/2.7K)- 4 . We take T C mb = 2.725. 

The redshif t of th e drag epoch Zd is well approximated by 
lEisenstein & Hul d 19981) 



Zd 



1291(fi m /i 2 



l + 0.659(fi m /i 2 ) - 828 



[l + bi(fi b ft 2 



(C9) 



where 

61 = 



0.313(fi m /i 
0.238(fi m /i 2 



[l + 0.607(fi m /i 2 ) ' 674 ] , (C10) 
(Cll) 



There are only four independent parameters among these five 
and n s is marginalized over in this study. Therefore, there are only 
three parameters left, {l a , R, z*}. CMB data are included in our 
analysis by adding the following term to the \ 2 of a given model 
with Apr = l a (z- t ) -302.35, Ap 2 = R(z*) - 1.728, and Ap 3 = 
z, - 1091.32: 

XCMB — 

where the inverse cov ariance matrix of (l a , R, z*) from WMAP7 
jKomatsu et alj l2010h is given by dWang. Chuang. & Mukherjeel 

Honl): 



Cov n \ rR = 



1.85710 
25.9289 
-1.14325 



25.9289 
5963.26 
-99.3185 



-1.14325 
-99.3185 
2.94429 



(C13) 
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